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摘 要: 采用 1956 一 2020 年 巩 勒 河上 游 昌 马 堡 水 文 站 逐 月 径流 量 资料 ,选取 基尼 系数 . 洛 伦 效 不 对 称 系数 等 8 个 


径流 年 内 分 配 特征 指标 ,结合 M-K 检 验 .RS 分 析 等 多 种 统计 方法 ,从 年 内 、 年 际 综合 分 析 玻 勒 河 上 游 径 流 演 变 规 
律 。 结 果 表 明 :(1) 玻 勒 河上 游 径 流量 年 内 分 配 呈 单 峰 型 分 布 ,主要 集中 在 汛 


, 除 绝对 变化 幅度 (AR ) 外 ,其 余年 


内 分 配 指标 均 呈 下 降 趋 势 ,表明 年 内 分 配 逐 渐 趋 向 均匀 ;(2) 年 径流 量 总 体 呈 增加 趋势 ,尤其 在 1997 年 发 生 突变 
后 ,倾向 率 加 大 ,年 径流 量 平均 值 较 突变 前 增加 59%;(3) 春 . 夏 .秋冬 四 季 径 流量 总 体 均 呈 上 升 趋势 ,其 中 夏季 倾 


向 率 最 大 ,冬季 倾向 率 最 小 ;20 世 纪 90 年 代 发 生 突变 后 ,径流 量 平 均值 较 突 变 前 均 明显 增 加 ,变化 幅度 排序 为 :秋季 
(76%)> 冬 季 (74% )> 夏 季 (589 )> 春 季 (45% ) 5 (4) 踊 勒 河上 游 年 和 四 季 径 流量 均 具 有 多 时 间 特 征 斥 度 周 期 , 且 第 一 


相似 内 陆 河 流域 径流 变化 规律 具有 指导 意义 。 


均 为 56 a; (5) 对 2022 一 2024 年 径流 量 预测 显示 ,下 勒 河上 游 径流 量 在 未 来 3 a 将 继续 呈 上 升 趋势 。 研 究 结 
果 为 准确 掌握 玻 勒 河 上 游 径 流 变 化 规律 及 特性 提供 了 科学 依据 ,对 玻 勒 河流 域 水 资源 的 可 持续 开发 利用 以 及 研究 


关键 词 : 径流 年 内 分 配 ; 年 际 变化 ; 径流 量 预 测 ;， 琉 勤 河 上 游 


径流 作为 水 系统 结构 的 重要 组 成 部 分 , 与 气候 
要 素 具 有 明显 的 相关 性 。 降 水 .气温 、 蔡 发 通过 
改变 陆 面 一 大 气 之 间 的 水 分 .能量 交换 及 传输 过 
程 ,进而 直接 改变 径流 变化 过 程 ,引起 水 资源 在 时 
间 和 空间 上 的 重新 分 配 ,加 剧 某 些 地 区 的 洪 洲 或 干 
旱灾 害 , 引 起 可 利用 水 资源 的 改变 ,进一步 对 生态 
系统 .自然 资源 .极端 气候 和 人 类 社会 产生 影响 。 
径流 变化 在 一 定 程度 上 会 导致 水 资源 量 与 水 生态 
系统 结构 发 生 改变 ,使 水 资源 合理 规划 与 高 效 管理 
面临 严峻 挑战 ,因此 ,全 面 分 析 变 化 环境 下 径流 演 
变 规 律 具 有 一 定 的 现实 意义 。 

为 厘清 变化 环境 下 径流 演变 特征 ,诸多 学 者 针 
对 一 些 典 型 流域 已 开展 了 大 量 工作 , 孙 甲 岚 等 汪 运 
用 M-K 检 验 .小 波 分 析 等 方法 系统 分 析 了 长 江 流域 
年 径流 量变 化 特征 ,Pan 等 5 基于 1951—2017 EX 
测 径流 数据 ,采用 累积 距 平和 M-K 检 验 对 黄河 流域 
径流 变化 进行 了 趋势 性 和 突变 性 分 析 , 邢 贞 相 等 9 
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采用 线性 倾向 .M-K 检 验 等 方法 分 析 了 挠 力 河 径流 
变化 特征 。 以 上 研究 大 多 集中 于 径流 量 的 年 际 变 
化 特征 分 析 ,但 随 着 气候 变 暖 和 人 类 活动 加 剧 , 径 
流 变 化 一 冰川 萎缩 一" 下游 河流 干 酒 等 一 系列 问 
题 ,使 得 径流 分 布 规律 发 生 显著 变化 。 径 流年 内 分 
布 情况 与 水 资源 合理 利用 、 区 域 防洪 抗旱 联系 紧 
密 , 部 分 学 者 针对 年 内 分 布 规 律 进行 了 相关 研究， 
涂 新 军 等 "通过 集中 度 .集中 期 AN E 
比 4 个 指标 ,探讨 东江 径流 的 年 内 变化 特征 ; 刘 永 婷 
等 “选取 基尼 系数 . 洛 伦 效 不 对 称 系数 2 个 年 内 分 
配 指标 ,定量 评价 淮河 上 游 径 流量 的 年 内 分 配 均匀 
性 。 总 的 来 说 ,径流 变化 主要 包括 数量 和 结构 2 种 
变化 "。 

玻 勒 河流 域 是 甘肃 省 河西 地 区 重要 水 源 地 之 
一 ,其 上 游 径流 量 的 变化 关系 着 中 .下游 地 区 水 资 
源 的 可 持续 性 利用 “”。 目 前 ,诸多 学 者 针对 玻 勤 
河上 游 径流 变化 规律 开展 了 大 量 工作 , 张 文 春 等 * 
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基于 昌 马 保水 文 站 的 年 平均 径流 量 , 分 析 了 玻 勒 河 
上 上 游 径流 量变 化 的 趋势 性 特征 ; 张 晓 晓 等 ”以 昌 马 
保水 文 站 的 实测 径流 量 为 研究 数据 ,基于 年 内 不 均 
SARL .变化 幅度 等 指标 ,探讨 了 玻 勒 河上 游 径流 
量 的 年 内 分 配 规 律 。 以 往 研 究 对 玻 款 河上 游 径流 
年 内 变化 .趋势 性 .突变 性 等 特征 分 析 较 为 全 面 ,但 
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夏季 .9 一 11 月 为 秋季 12 AKA RB, 
1.3 研究 方法 

1.3.1 年 内 分 配 分 析 方法 为 了 定量 分 析 玻 勒 河 上 
游 径流 演变 规律 ,在 基于 月 (或 季 ) 占 年 径流 的 百分比 
等 方法 分 析 的 基础 上 ,选取 年 内 不 均匀 系数 (C,) 2 、 
SCA WD BUC.) RRE) RPD) 


在 年 内 分 配 指标 选取 、 特 征 分 析 方 法 方面 还 可 进 一 
步 补充 完善 。 因 此 ,本 文 从 年 内 、 年 际 两 方面 出 发 ， 
采用 M-K 检 验 法 .RS 分 析 法 .小波 分 析 等 统计 学 方 
法 综合 分 析 玻 勒 河上 游 径流 年 内 分 配 情况 .趋势 
性 ` 变 异性、 持续 性 及 周期 性 ,并 运用 灰色 预测 模 
型 .R/S 灰色 预测 模型 预测 未 来 年 径流 量 ,系统 研究 
了 径流 时 间 变 化 规律 ,为 进一步 明晰 环境 变化 下 内 
陆 河 流域 径流 变化 影响 因素 黄 定 基础 。 


1 数据 与 方法 


1.1 研究 区 概况 

芒 勒 河流 域 是 甘肃 省 三 大 内 陆 河 流域 之 一 ， 
以 昌 马 水 库 、 双 塔 保水 库 作 为 上 .中 、 下 游 的 分 界 
点 057 。 跤 勒 河上 游 河 段 全 长 328 km ,流域 面积 约 
为 10961 km ,属于 干 冷 .多 风 的 高 原 大 陆 性 气候 ,多 
年 平均 径流 量 为 10.28x10 mi, 冰 川 、 积 雪 和 冻 土 分 
布 面积 较 广 趾 ,区 域内 设 有 昌 马 保水 文 站 (图 1)。 
1.2 数据 来 源 

数据 来 自 于 甘肃 省 酒泉 水 文 站 玉 勒 河上 游 的 控 


制 水 文 站 一 一 昌 马 保水 文 站 1956 一 2020 年 的 实测 

月 径流 数据 。 文 中 将 3 一 5 月 划分 为 春季 、6 一 8 月 为 
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图 1 研究 区 位 置 示意 图 
Fig.1 Schematic diagram of study area 


洛 伦 兹 不 对 称 系数 (S)' ”基尼 系数 (GI)'Y、 相 对 变 
化 幅度 (C.) 2 和 绝对 变化 幅度 (AR) 2 8 个 径流 年 内 
分 配 特征 指标 ,从 不 同 角 度 分 析 径 流年 内 变化 
1.3.2 年 际 变 化 分 析 方 法 ”运用 R/S 分 析 法 、M-K 检 
验 、 滑 动 : 检 验 、 线 性 倾向 佑 计 和 小 波 分 析 等 方法 ， 
从 年 际 方面 分 析 玻 勒 河上 游 径流 量 持续 性 、 突 变 
性 .趋势 性 以 及 周期 性 特征 ?? 。 

1.3.3 径流 预测 分 析 方 法 ”灰色 预测 模型 是 一 种 优 
于 线性 回归 的 单 序列 预测 模型 ,其 序列 表达 式 为 一 
阶 线性 微分 方程 于。 由 于 径流 变化 具有 突变 性 ， 
为 了 提高 预测 精度 ,提出 了 RS 灰色 预测 模型 于 ?9 ， 
其 计算 原理 为 :(1) 根据 R/S 分 析 法 ,确定 序列 平均 
循环 周期 1; (2) 以 N=t-1 为 研究 周期 建立 灰色 预测 
模型 ,预测 N+1 年 径流 量 ;(3) 对 模型 进行 残 差 检 
验 , 若 通过 , 则 依次 类 推 ,继续 预测 下 一 年 径流 量 。 
ARIMA 模 型 是 将 ARMA 模型 与 差分 运算 结合 的 组 合 
预测 模型 ,又 称 为 差分 整合 移动 平均 自 回归 模型 ” 
iA ARIMA (p,d.q) ,其 中 pb 为 自 回归 项 数 ,d 为 序列 
平稳 时 的 差分 次 数 ,q 为 滑动 平均 项 数 ™。 


2 结果 与 分 析 


2.1 径流 年 内 分 配 特征 

2.1.1 年 内 分 配 情况 ” 芷 勒 河流 域 是 典型 的 内 陆 河 
流域 ,其 主要 补给 来 源 为 降水 和 冰川 融雪 。 由 图 2a 
可 知 , 昌 马 堡 站 径流 量 呈 单 峰 型 分 布 ,冬季 处 于 低 
值 ;春季 气温 回升 ,降水 和 融雪 增加 ,径流 量 呈 上 升 
趋势 ;夏季 汛期 来 临 , 径 流量 急剧 上 升 , 在 8 月 达到 
最 大 值 ; 秋 季 和 气温 降低 ,径流 量 呈 缓慢 下 降 趋 势 , 至 
12 月 回落 到 常态 。 结 合 图 2b 可 知 ,在 各 年 代 内 分 配 
上 ,径流 主要 集中 在 汛期 7 一 8 月 , 占 多 年 平均 值 的 
49.23% ~53.29% ; 2010s 、2000s 、1980s .1950s 的 月 最 
大 径流 主要 出 现在 7 月 ,分 别 占 多 年 平均 值 的 
24.20% , 24.43% , 23.74% All 27.97% ; 1990s , 1970s , 
1960s 的 月 最 大 径流 主要 出 现在 8 月 ,分 别 占 多 年 平 
均值 的 25.35% .26.60% 和 24.74%。 
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图 2 玻 勒 河上 游 径 流年 内 分 配 过 程 


Fig. 2 Annual distribution process ofrunoff in the upper reaches of the Shule River 


由 表 1 可 知 , 琉 勒 河 上 游 多 年 平均 径流 量 ， 
10.28x10° m ,各 时 段 径 流 呈 持续 增加 状态 ,2010s 时 
段 径 流量 较 1950s 增 加 了 1 倍 多 。 全 年 均 保 持 上 升 
趋势 ,其 中 ,四 季 占 各 时 段 均值 的 具体 排序 为 :夏季 
(56.50% ~62.91% ) > 秋季 (17.60% ~20.39% ) > 春季 
(11.41%~15.48% )> 冬 季 (7.25%~8.97% ) 。 总 体 上 ， 
基于 径流 上 升 的 大 趋势 ,全 年 径流 集中 在 夏秋 季 ， 


\ 


日 一 直 存 在 夏季 径流 比例 偏 大 、 春 季 径 流 比 例 偏 
的 现象 。 

2.1.2 年 内 分 配 特 征 受气 候 变 化 和 人 类 活动 的 影 
响 , 琉 勒 河上 游 水 循环 发 生 了 明显 的 季节 性 变化 ， 
进而 导致 流域 径流 年 内 分 配 不 均匀 性 发 生 改 变 。 
由 表 2 可 知 ,C,、 GCC、D 分 别 以 - 0.015.(10aj 
一 0.006(10a)'、 一 0.006(10a)"'、 一 1.053.(10a)"' 的 速 


= 


表 1 径流 年 内 各 季 分 配 


Tab. 1 Seasonal distribution of runoff during the year 


een 春季 (3 一 5 月 ) 夏季 (6 一 8 月 ) 秋季 (9 一 11 月 ) 冬季 (12 月 一 次 年 2 月 ) 
时 自 fa 径流 量 百分比 ”径流 量 百分比 ”径流 量 。 百分比 Aie ”百分比 
/10°m /% /10°m* 1% /10°m /% /10°m* 1% 
1950s 7.37 0.84 11.41 4.63 62.91 1.36 18.43 0.53 7.25 
1960s 8.03 1.24 15.48 4.79 59.60 1.41 17.60 0.59 7.32 
1970s 8.92 1.27 14.22 5.23 58.61 1.71 19.14 0.72 8.04 
1980s 9.22 1.38 15.00 5.41 58.68 1.66 18.02 0.76 8.29 
1990s 8.92 1.23 13.83 5.27 59.08 1.68 18.83 0.74 8.26 
2000s 12.72 1.80 14.14 7.19 56.50 2.59 20.39 1.14 8.97 
2010s 14.87 1.78 11.94 8.84 59.47 3.00 20.18 1.25 8.41 
多 年 平均 值 10.28 1.42 13.82 6.03 58.64 1.98 19.26 0.85 8.28 
R2 朴 勒 河上 游 径流 年 内 分 配 特征 指标 
Tab.2 Intra-year distribution characteristics of runoff in the upper reaches of the Shule River 
统计 值 G; G, G D GI $ Ci AR 

倾向 率 /10a)” -0.015 -0.006 -0.006 -1.053 -0.009 -0.006 -0.643 0.255 

最 大 值 1.27 0.66 0.66 276.48 0.49 1.09 32.20 5.66 

最 大 值 出 现年 份 1976 1976 1958 1987 1991 1975 1958 2016 

最 小 值 0.68 0.26 0.39 181.23 0.29 0.72 6.81 0.70 

最 小 值 出 现年 份 1973 1984 2003 1969 2020 2019 1973 1956 

平均 值 0.93 0.36 0.52 197.42 0.39 0.84 12.16 2.54 
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率 下 降 , 在 不 考虑 2016 年 和 1956 年 极 值 的 情况 下 ， 
Ca MAR 也 呈 下 降 趋 势 , 即 年 内 最 大 月 径流 量 和 最 
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变 ;综合 考虑 以 上 2 种 突变 分 析 方 法 检验 结果 ,可 认 
为 疏 勒 河 上 游 年 径流 量 在 1997 年 发 生 突变 。 由 表 3 


小 月 径流 量 的 差 值 不 断 减 小 ,从 不 均匀 性 、 集 中 度 、 
变化 幅度 等 方面 来 看 , 玻 勒 河上 游 年 内 分 配 趋向 均 
名;65 a 中 GI 均 分 布 在 0.29~0.49 之 间 , 其 中 有 40a 
GI 值 小 于 0.4, 自 2010 年 以 来 ,GI 值 一 直 小 于 0.4, 表 
明 玻 勒 河上 游 年 内 分 配 趋 向 均匀 ;65 a 的 S 大 都 小 
于 1, 表 明 径 流年 内 分 配 不 均匀 的 原因 可 能 是 径流 
量 较 小 的 月 份 造成 的 。 
2.2 径流 持续 性 .突变 性 分 析 

Hurst 指数 可 以 定量 表征 时 间 序 列 的 长 期 相关 
性 ,其 值 大 小 可 作为 时 间 序 列 是 否 存在 趋势 性 成 分 
的 依据 ,一 般 通 过 R/S 分 析 计 算 Hurst 指 数值 。 结 合 
图 3 和 表 3 可 知 , 玻 勒 河 上 游 年 径流 量 和 四 季 径 流 
量 的 Hurst 指数 均 大 于 0.5 ,表现 为 正 持续 性 ,其 中 ， 
除 夏 季 径 流量 的 Hurst 指数 为 0.7353 ,年 径流 量 A 
季 径 流量 .秋季 径流 量 和 冬季 径流 量 的 Hurst 指数 
均 远 远大 于 0.5 ,表现 出 很 强 的 正 持 续 性 ,表明 未 来 
臣 勒 河上 游 径流 量 将 继续 保持 上 升 趋势 。 

H M-K 4 g (| 4a) PJ EN , UF FI UB FE 95% E 
信 区 间 内 存在 交点 (1997 年 ) ,表明 1997 EA a Tay 
上 游 年 径流 量 的 突变 点 ;采用 10 a 滑 动 : 检 验 (图 
4b) 分 析 , 得 到 年 径流 量 在 20 世 纪 90 年 代 末 发 生 突 


3.5 F (a) 


3.0 F y=0.82793x-0.80576 / 
R°=0.82519 


25 F 
2.0 F 


In (R/S) 


IS F 
1.0 上 
05 f 


0.0 上 


AY, RD EREE 夏季 秋季 和 冬季 的 突变 年 
份 分 别 为 1999 年 .1998 年 .1997 年 和 1999 年 。 
2.3 径流 趋势 性 分 析 

由 图 5a 可知 , 玻 勒 河上 游 年 平均 径流 量 ， 
10.28 10° m°, 最 大 值 为 17.44x10xsm ,出 现在 2017 
年 ,最 小 值 为 4.13x10sm ,出 现在 1956 年 。 结 合 表 4 
进一步 分 析 ,年 径流 量变 化 趋势 通过 了 0.05 显著 性 
水 平 ,总 体 上 以 1.21x10sm.(10a 六 的 速率 显著 上 升 ， 
且 年 径流 量 在 突变 后 较 突 变 前 增加 了 59% 。 由 图 
5b 可知 ,1997 年 前 ,年 平均 径流 量 累 积 距 平 值 呈 下 
降 趋 势 ,年 平均 径流 量 呈 下 降 趋 势 ;1997 年 后 ,年 平 
均 径 流量 累积 距 平 值 呈 上升 趋势 ,年 平均 径流 量 呈 
逐步 上 升 趋势 日 上 升 速率 较 快 ,这 与 突变 检验 结 
相 一 致 ,再 一 次 印证 玖 勒 河上 游 年 径流 量 在 1997 年 
发 生 了 突变 。 

结合 图 6 和 表 4 分 析 可 知 , 琉 勒 河 上 游 四 季 丛 
流量 均 通 过 了 0.05 的 显著 性 检验 ,总 体 均 呈 现 上 升 
趋势 。 其 中 ,夏季 变化 幅度 最 大 ,倾向 率 为 0.66， 
(10a) ,秋季 次 之 ,倾向 率 为 0.29. (10a)' ,春季 再 次 ， 
倾向 率 为 0.14. 10a) " ,冬季 增幅 最 小 ,气候 倾向 率 
为 0.13* 40a)' 。 由 表 4 可 知 , 夏 季 、 秋 季冬 季 突 变 


3:0 F (b) 
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N 
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F 径 流量 的 R/S 分 析 结 果 


Fig. 3 R/S analysis results of annual runoff in the upper reaches of the Shule River 


3 HMA LHS 


径流 量 持续 性 突变 性 分 析 结 果 


Tab.3 Analysis results of annual and four-season runoff persistence and mutation in the upper reaches of the Shule River 


1956—2020 年 径流 序列 年 
持续 性 分 析 Hurst 指 数 0.8279 
WER 0.8252 
突变 性 分 析 突变 年 份 1997 


春季 夏季 秋季 冬季 
0.8127 0.7353 0.8155 0.9422 
0.9690 0.6802 0.8846 0.9738 

1999 1998 1997 1999 
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Fig.4 Analysis on the mutation of annual runoff in the upper reaches of Shule River 


(a) 


y=0.12077x-229.8104 
R°=0.51675 


年 径流 量 /108ms 
= 


—— 年 平均 径流 量 --- 突变 前 径流 量 
… 突变 后 径流 量 -,- 线性 (年 平均 径流 量 ) 
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年 份 


图 5 世 勒 河上 游 年 径流 量 趋势 线 和 距 平 


#E3E/108 m3 


累积 距 平 /108m 
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年 份 


Fig. 5 Trend lines and anomalies of annual runoff in the upper reaches of the Shule River 


后 径流 量 的 变 差 系数 均 小 于 突变 前 ,表明 突变 后 径 
流量 变异 性 较 小 , 且 四 季 突 变 后 的 径流 量 均 值 较 突 
变 前 大 幅 增 加 ,具体 排序 为 :秋季 (76% ) > 冬季 
Di 2 

4 径流 周期 性 分 析 

由 图 7a 方 差 图 可 知 ,下 勒 河上 游 年 径流 量 存在 
3 个 主 周期 ,分 别 为 56a.28 a 和 10 a, 其 中 56 a 为 第 
一 主 周期 。 结 合 图 7b 实 部 图 进一步 分 析 , 对 于 56 a 
的 特征 时 间 尺 度 , 存 在 2 个 明显 的 丰 水 期 中 心 和 1 
个 村 水 期 中 心 ,对 应 年 份 分 别 为 1974 年 .2010 年 和 
1989 年 ,经 历 了 1 个 丰 枯 转化 期 。28 a 尺度 上 存在 3 
个 丰 水 期 中 心 和 3 个 枯水期 中 心 ,经 历 了 3 个 丰 枯 
转化 期 ,其 中 丰 水 期 中 心 为 1967 年 .1985$ 年 和 2004 
年 ,枯水期 中 心 为 1975 年 .1994 年 和 2013 年 。10 a 
尺度 上 经 历 了 大 约 10 个 丰 枯 转化 期 , 丰 枯 交替 频 
繁 。 由 表 5 TA, Gi LEY ES TN MA 


尺度 特征 ,春季 的 第 一 二 和 第 三 主 周期 分 别 为 56 
a28 a 和 20 a, 夏 季 的 第 一 、 第 二 和 第 三 主 周期 分 别 
为 56a、28a 和 9 a, 秋 季 径 流 的 第 一 、 二 主 周期 分 别 
为 56a 和 13 a, 冬 季 径 流 的 第 一 主 周期 为 56 a。 
25 径流 量 预测 分 析 

用 1956 一 2021 年 玩 勒 河上 游 实测 年 径流 量 资 
料 ,分 别 构建 灰色 预测 模型 .R/S 灰色 预测 模型 、 
ARIMA 模型 ,以 2019 一 2021 年 为 验证 阶段 ,预测 
2022—2024 年 径流 量 ,分 析 结 果 见 表 6。 综 合 考 虑 3 
种 模型 的 预测 效果 ,最 终 以 综合 预测 值 作 为 2022 一 
2024 年 径流 量 预测 值 ,分 别 为 15.49x10:m’、15.53x 
10° m?,15.67X10°m’, 


3 讨论 


BTA EWE EF AER , 受 人 类 活动 干扰 较 
小 ,径流 量 总 体 呈 上 升 趋势 ,本 研究 基于 1956 一 
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图 6 下 勒 河上 游 四 季 径 流量 趋势 分 析 
Fig.6 Analysis on the trend of runoff in four seasons in the upper reaches of Shule River 
4 琉 勒 河上 游 径流 量 趋势 分 析 
Tab.4 Analysis of runoff trend in the upper reaches of the Shule River 
统计 特征 趋势 检验 突变 点 前 突变 点 后 
1956 一 2020 年 Oi — a rE ih vE 
径流 序列 os 项 向 率 M-K 显著 | 人 me 
aaa /10°m? 系数 倾向 率 。 MA 检验 值 。 显著 性 /10°m* 系数 /10°m* 系数 
年 10.28 0.31 1.21 6.0633 是 8.50 0.21 13.54 0.18 
春季 1.42 0.26 0.14 5.8482 是 1.24 0.18 1.80 0.20 
夏季 6.03 0.34 0.66 5.0386 是 5.04 0.27 7.96 0.22 
秋季 1.98 0.40 0.29 5.7463 是 1.56 0.28 2.75 0.25 
冬季 0.85 0.33 0.13 7.6825 是 0.69 0.18 1.20 0.14 
(a) 方差 图 (b) 小 波 实 部 图 SHEE 
1000 n 
4 
800 3 
= 1 
ic 
人 3 e 
A 400 Ke -2 
200 2 K 
5 
0 
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Fig. 7 Wavelet period analysis of annual runoff in the upper reaches of the Shule River 
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Tab. 5 Annual and four-season cycles in the upper 


reaches of the Shule River 


1956—2020 年 径流 序列 第 一 主 周期 第 二 主 周期 第 三 主 周期 


年 56a 28a 10a 
春季 56a Ba 20 a 
夏季 56a Ba 9a 
秋季 56a 13a - 
冬季 56a 一 一 


2020 年 昌 马 保水 文 站 逐 月 径流 量 分 析 径流 演变 特 
征 ,发 现 研 究 时 段 内 , 春 、 夏 .秋冬 四 季 年 径流 量 总 
体 均 呈 上 升 趋势 ,其 中 夏季 上 升 趋势 最 高 ,冬季 上 
升 趋势 最 低 。 丁 宏伟 等 号 以 傅 里 时 公式 为 基础 的 
频谱 分 析 方 法 发 现 , 玻 勒 河上 游 1952 一 1999 年 径流 
量 以 6a 11a 和 3a 周 期 变化 最 为 显著 ;本 研究 通过 
小 波 分 析 方 法 发 现 ,1956 一 2020 年 径流 量 最 显著 的 
周期 为 56 a、28 a 和 10a。 由 于 选取 研究 时 段 的 序 
列 长 度 不 一 致 或 采用 方法 不 同 ,会 导致 分 析 结 果 有 
所 不 同 ; 且 径流 具有 随机 性 、 突 变性. 非 线 性 等 复杂 


表 6 GRA Lee 


特性 外 , 仅 用 一 种 方法 分 析 会 造成 结果 的 局 限 性 和 
不 确定 性 ,不 利于 全 面 认识 环境 变化 下 径流 演变 
规律 。 
此 外 ,基于 朴 勒 河上 游 径 流 的 年 内 不 均匀 系数 
C) 完全 调节 系数 (C) EPEC) 、 集 中 期 (D)、 
相对 变化 幅度 (C,) 和 绝对 变化 幅度 (AR)6 个 指标 ， 
增加 了 洛 伦 效 不 对 称 系数 ($S) 和 基尼 系数 (GI)2 个 
虽 标 分 析 径 流量 对 均匀 性 的 贡献 程度 后 , 发现 除 绝 
对 变化 幅度 外 ,其 余 分 配 指标 均 呈 下 降 趋 势 。 在 不 
考虑 极 值 时 ,绝对 变化 幅度 也 呈 下 降 趋 势 ,进一步 
说 明年 内 分 配 过 程 逐渐 趋向 均匀 。 采 用 灰色 预测 
模型 .R/S 灰色 预测 和 ARIMA 模型 分 别 预 测 未 来 玻 
勒 河上 游 年 径流 量 ,发 现 灰色 预测 模型 的 相对 误差 
最 小 ,但 若 遇 到 突变 年 份 , 就 会 造成 较 大 的 误差 。 
考虑 径流 变化 具有 周期 性 ,通过 R/S 分 析 方 法 确定 
循环 周期 7, 以 7-1 为 研究 序列 进行 灰色 预测 ,但 是 
R/S 灰色 预测 模型 与 研究 序列 的 波动 性 和 跳跃 性 有 
关 , 跳 跃 性 越 强 、 波 动 性 越 大 ,其 预测 精度 就 越 小 。 


径流 量 预测 结果 


Tab.6 Prediction results of annual runoff in the upper reaches of the Shule River 


灰色 预测 模型 R/S 灰色 预测 模型 ARIMA 模型 

qp FNA tr = n o - 2 预测 值 i ae ee 

/10°m° i : à 108m’ 误差 /% 

/10°m? 1% /10° m? 1% /10°m? 1% 

2019 16.91 13.92 -17.69 16.28 -3.69 15.73 -6.96 15.31 -9.46 
2020 12.25 14.04 14.59 16.81 37.23 15.88 29.61 15.58 27.18 
2021 14.03 14.16 0.89 17.06 21.57 16.02 14.16 15.75 12.26 
2022 14.33 16.89 15.26 15.49 
2023 14.45 16.74 15.39 15.53 
2024 14.57 16.94 15.51 15.67 
平均 相对 误差 /% 11.06 20.83 16.91 16.30 


因此 ,本 文采 用 ARIMA 等 单 变量 预测 模型 ,主要 是 
根据 径流 自身 变化 规律 进行 预测 分 析 ,未 能 考虑 降 
雨 .气温 等 气象 条 件 对 径流 变化 的 影响 ,使 得 预测 
值 与 实测 值 之 间 存 在 一 定 的 误差 , 故 今后 需要 继续 
开展 多 变量 预测 模型 等 相关 内 容 的 研究 。 


4 结论 

本 研究 基于 巩 勒 河 上 游 昌 马 堡 水 文 站 1956 一 
2020 年 月 径流 数据 ,研究 径流 变化 特征 ,分析 结果 
如 下 : 


(1) BP EWE TALE AL IR OBE MT 
称 系数 等 年 内 分 配 指标 均 呈 下 降 趋 势 , 表 明年 内 分 
配 过 程 逐渐 趋向 均匀 , 主要 是 因为 气温 升 高 导致 枯 
水 季节 径流 量 增加 ,进而 在 一 定 程 度 上 调节 了 径流 
年 内 分 配 均匀 度 。 

(2) 年. 春 . 夏 . 秋 、 冬 径流 量 在 20 世 纪 90 年 代 
末 发 生 突变 后 ,径流 量 平均 值 较 突变 前 分 别 增加 了 
59% 、45% 、58% .76% 、74% ,这 与 近年 来 降水 增加 、 
气温 升 高 造成 冰雪 径流 量 增加 密切 相关 。 

(3) 玻 勒 河 上 游 年 和 四 季 径 流量 均 具 有 多 时 间 
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特征 尺度 周期 。 其 中 ,第 一 主 周期 均 为 56 a, 年 径流 
量 .春季 径流 量 和 夏季 径流 量 的 第 二 主 周期 均 为 
28 aco 

(4) BHA oe LEE E 四季 径 流 
量 均 呈 上 升 趋势 ,经 综合 预测 , A eT E WE it E 
在 未 来 3 a 继 续 呈 上 升 趋势 ,预测 值 分 别 为 15.49x 
10° m’ 15.53x10° m’ 15.67X10°m’, 
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Analysis of runoff variation and forecast in the upper reaches of the Shule River 


JIA Ling, ZHANG Baizu’, NIU Zuirong', SUN Dongyuan ， 
SUN Kaiyue, WNAG Lujun' 
(1. College of Water Resources and Hydropower Engineering, Gansu Agricultural University, Lanzhou 730070, 
Gansu, China; 2. Jiuquan Hydrological Station of Gansu Province, Jiayuguan 735100, Gansu, China) 


Abstract: A quantitative analysis of runoff variation characteristics is the basis for finding out its impact 
mechanism, and it is also an important basis for accurately evaluating water resources. This study used eight 
annual distribution indexes, such as the Gini coefficient and Lorentz asymmetry coefficient, combined with the 
M-K test, R/S analysis, and other mathematical statistics methods, to comprehensively analyze the process of the 
runoff of the upper reaches of the Shule River from the annual and interannual levels. The results showed that the 
following: (1) during the year, a unimodal distribution was shown, which was mainly concentrated in the flood 
season concentration degree (D), inhomogeneous coefficient (Cy), Gini coefficient (GI), relative change range 
(Cm) decreased respectively at the rate of- 1.053 - (10a) ', - 0.015 - (10a) ', - 0.009 - (10a) ',- 0.643 - (10a) ', 
complete adjustment coefficient (C;), Lorentz asymmetry coefficient (S), concentration period (Cj all decreased 
at the rate of -0.006 -(10a)'. The distribution tended to be uniform gradually according to the changing trend of 
its characteristic index during the year. (2) The annual runoff showed an overall increasing trend, but an abrupt 
change occurred in 1997, and the average annual runoff afterward increased by 59% compared before the change. 
(3) The runoff generally showed an upward trend in all four seasons, among which the trend rate was the greatest 
in the summer and the smallest in the winter. After the abrupt change in the 1990s, the average runoff increased 
significantly compared to before the abrupt change, and the change range was in the following order: autumn 
(76%) > winter (74%) > summer (58%) > spring (45%). (4) The annual and seasonal runoff in the upper reaches 
of the Shule River had multiple time characteristic scale periods, and the first major cycle were both 56 years. (5) 
According to the predicted runoff from 2022 to 2024, it was found that the runoff of the upper reaches of the 
Shule River will continue to increase in the next three years. These results provided a scientific basis for 
accurately grasping the runoff variation laws and characteristics of the upper reaches of the Shule River and had 
certain guiding significance for the sustainable development and utilization of water resources in the Shule River 
Basin and the study of the runoff variation laws in inland river basins. 


Keywords: runoff seasonal distribution; inter-annual variation; runoff forecast; upper reaches of the Shule River 


